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DRIFT-FREE VIDEO ENCODING AND DECODING METHOD , AND CORRESPONDING DEVICES 



FIELD OF THE INVENTION 

The present invention relates to an encoding method for the compression of an 
original video sequence divided into successive groups of frames (GOFs) and to a 
corresponding decoding method. It also relates to corresponding encoding and decoding 
devices. 

BACKGROUND OF THE INVENTION 



enabled new applications and services. Many of them not only require coding efficiency but 
also enhanced functionality and flexibility in order to adapt to varying network conditions 
and terminal capabilities. Scalability answers these needs. Current video compression 
standards often use so-called hybrid solutions, based on a predictive scheme where each 
frame is temporally predicted from a reference frame (the prediction options being : zero 
value prediction, for the intra frames or I frames, forward prediction, for the P frames, or bi- 
directional prediction, for the B frames) and the obtained prediction error is spatially 
transformed to get advantage of spatial redundancies. From MPEG-2 to MPEG-4, standard- 
based scalable solutions have then been proposed. They rely on the generation of a base 
layer, containing the lowest spatial, temporal and/or SNR resolution version of the original 
video sequence, and one or several enhancement layers allowing (if transmitted and decoded) 
a spatially, temporally and/or SNR refined reconstruction. A short-coming of these layer- 
based scalability schemes comes however from their lack of coding efficiency. 



dimensional (3D) subband coding, which are able to generate embedded bitstreams. Thanks 
to their multi-resolution analysis structure, scalability is inherent to these schemes and does 
not weaken their intrinsic coding efficiency. In a 3D subband codec such as described for 
example in "A fully scalable 3D subband video codec", "Proceedings of the International 
Conference on Image Processing (ICIP2001), vol.2, 2001, pp.1017-1020, the embedded 
bitstream is fully scalable and can be decoded at any spatial and temporal resolutions, and 
with any desired SNR quality, simply by truncation at known locations. In such a scheme, 



The growth of the Internet and advances in multimedia technologies have 



A different approach has been proposed with techniques such as three- 
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successive groups of frames (GOFs) are processed as a 3D structures and spatio-temporally 
filtered in order to compact the energy in the low frequencies, a motion compensation being 
also provided in order to improve the overall coding efficiency. The 3D subband structure is 
depicted in Fig.l : the illustrated 3D wavelet decomposition with motion compensation is 
5 applied to a group of frames (GOF), and this current GOF is first motion-compensated (MC), 
in order to process sequences with large motion, and then temporally filtered (TF) using Haar 
wavelets (the dotted arrows correspond to a high-pass temporal filtering, while the other ones 
correspond to a low-pass temporal filtering). After the motion compensation operation and 
the temporal filtering operation, each temporal subband is spatially decomposed into a 

1 0 spatio-temporal subband, which finally leads to a 3D wavelet representation of the original 
GOF, three stages of decomposition being shown in the example of Fig. 1 (L and H = first 
stage ; LL and LH = second stage ; LLL and LLH = third stage).The well known SP1HT 
algorithm, extended from 2D to 3D, is chosen in order to efficiently encode the final 
coefficient bit-planes with respect to the spatio-temporal decomposition structure. 

15 As it is implemented now, a 3D subband codec applies the 

motion-compensated (MC) spatio-temporal analysis at the full original resolution at the 
encoder side. Spatial scalability is achieved by getting rid of the highest spatial subbands of 
the decomposition. However, when motion compensation is used in the 3D analysis scheme, 
this method does not allow a perfect reconstruction of the video sequence at lower resolution, 

20 even at very high bit-rates : this phenomena, referred to as drift in the following description, 
lowers the visual quality of the scalable solution compared to a direct encoding at the 
targeted final display size. As explained in the document "Multiscale video compression 
using wavelet transform and motion compensation", P.Y.Cheng and al., Proceedings of the 
International Conference on Image Processing (ICIP95), Vol.1, 1995, pp.606-609, this drift 

25 comes from the order of wavelet transform and motion compensation that is not 

interchangeable. Indeed, when a frame (A) is synthesized at a lower resolution (a), the 
following operation is applied : 

a = DWT L (L)+ MC[DWT L (H)] 

= DWT L (A)+ [MC[DWT L (H)] - DWT L (MC[H])] (1) 
30 where DWT L denotes the resolution downsample using the same wavelet filters as in the 3D 
analysis. In a perfect scalable solution, one wants to have: 

a = DWTl (A) (2) 
The remaining part of the expression (1) therefore corresponds to the drift. It can be noticed 
that, if no MC is applied, the drift is removed. The same phenomena happens (except at the 
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image borders) if a unique motion vector is applied to the frame. Yet, it is known that MC is 
unavoidable to achieve a good coding efficiency, and the likelihood of a unique global 
motion is small enough to eliminate this particular case in the following paragraphs. 

Some authors, such as J.W.Woods and al in the document "A resolution and 

5 frame-rate scalable subband/wavelet video coder", IEEE Transactions on Circuits and 
Systems for Video Technology, vol.1, n°9, September 2001, pp.1035-1044, get rid of this 
drift to achieve good spatial scalability by different means. However, in said document, the 
described scheme, in addition to being quite complex, implies the sending of an extra 
information (the drift correction necessary to correctly synthesize the upper resolution) in the 

10 bitstream, thus wasting some bits (the solution described in the document <( Multiscale video 
compression. . ." avoids this bottleneck but works on a predictive scheme and is not 
transposable to the 3D subband codec). 



SUMMARY OF THE INVENTION 
15 It is therefore an object of the invention to propose a solution avoiding these 

drawbacks. 

To this end, the invention relates to a video encoding method for the 
compression of an original video sequence divided into successive groups of frames (GOFs), 
said method comprising the steps of : 
20 (1) generating from the original video sequence, by means of a wavelet 

decomposition, a low resolution sequence including successive low resolution GOFs ; 
(2) performing on said low resolution sequence a low resolution decomposition, 

by means of a motion compensated spatio-temporal analysis of each 
low resolution GOF ; 

25 (3) generating from said low resolution decomposition a full resolution sequence, 

by means of an anchoring of the high frequency spatial subbands resulting from the wavelet 
decomposition to said low resolution decomposition ; 

(4) coding said full resolution sequence and the motion vectors generated during 

the motion compensated spatio-temporal analysis, for generating an output coded bitstream. 
30 The proposed solution is remarkable in the sense that the global structure of 

the decomposition tree in the 3DS analysis is preserved and no extra information is sent to 
correct the drift effect (only the decomposition/reconstruction mechanism is changed). If no 
motion estimation/compensation is performed at full resolution, it is a low-cost solution in 
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terms of complexity. If motion compensation is introduced in the high spatial subbands, a 
better coding efficiency is provided. 

The invention also relates to a corresponding decoding method, comprising the 

steps of : 

5 (1) decoding said input coded bitstream for generating a decoded full resolution 

sequence and associated decoded motion vectors ; 

(2) in said decoded full resolution sequence, separating the decoded high 
frequency spatial subbands and the decoded low resolution decomposition ; 

(3) generating from said decoded low resolution decomposition, by means of a 
10 motion compensated spatio-temporal synthesis, a decoded low resolution 

sequence ; 

(4) reconstructing from said decoded low resolution sequence and the decoded 
high frequency spatial subbands an output full resolution sequence corresponding to the 
original video sequence. 

15 The invention also relates to an encoding device and a decoding device 

provided for implementing said encoding method and said decoding method respectively. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The invention will now be described in a more detailed manner, with reference 
20 to the accompanying drawings in which : 

Fig.l shows a 3D subband decomposition ; 

Fig.2 illustrates a motion-compensated temporal analysis at the lowest 

resolution ; 

Fig.3 depicts an embodiment of an encoding scheme according to the 

25 invention ; 

Fig.4 depicts an embodiment of a decoding scheme corresponding to the 
encoding scheme of Fig.3 ; 

Fig. 5 illustrates the reordering of the high spatial subbands (for a forward 
motion compensation) ; 
30 Fig.6 depicts another embodiment of an encoding scheme according to 

the invention. 



DETAILED DESCRIPTION OF THE INVENTION 
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The proposed solution (i.e. a spatial scalability with no drift in a motion 
compensated 3D subband codec) is now explained with reference to its two main steps : (a) 
motion compensation at the lowest resolution, (b) encoding the high spatial subbands. 

First in order to avoid drift at lower resolutions, Motion Compensation (MC) 

5 is applied at this level. Consequently, as illustrated in Fig.2, one first downsizes (reference d) 
the GOF using wavelet filters, and the usual 3D subband MC-decomposition scheme is then 
applied to this downsized GOF instead of the full-size GOF. In Fig.2, the temporal subbands 
(L 0 ,d, H 0)d ) and (Li, dj Hi^j) are determined according to the well-known lifting scheme (H is 
first defined from A and B, and then L from A and H), and the dotted arrows correspond to 

10 the high-pass temporal filtering, the continuous ones to the low-pass temporal filtering, and 
the curved ones (between low frequency spatial subbands A of the frames of the sequence, 
referenced A 0 ,d, Ai >d , A 2 ,d, A3,d, or between low frequency temporal subbands L, referenced 
L 0> a and Li td ) to the motion compensation (it may be noticed that a side effect of this method 
is the reduction of the amount of motion vectors to be sent in the bitstream, which saves up 

15 some bits for texture coding). Before transmitting the subbands to a tree-based entropy coder 
(for instance to a 3D-SPIHT encoder such as described for instance in the document "Low 
bit-rate scalable video coding with 3D set partitioning in hierarchical trees (3D-SPIHTy\ B.J. 
Kim and al. IEEE Transactions on Circuits and Systems for Video Technology, vol.10, n°8, 
December 2000, pp. 1374-1387), one puts the high spatial subbands that allow the 

20 reconstruction of the full resolution. The final tree structure looks very similar to that of a 3D 
subband codec such as the one described in the document "A fully scalable 3D subband 
video codec", IEEE Conference on Image Processing (ICIP2001), vol.2, pp.1017-1020, 
Thessalonild, Greece, October 7-10, 2001, and so a tree-based entropy coder can be applied 
on it without any restriction, as described in the new encoding scheme of Fig.3, where the 

25 references are the following (for a frame of the full resolution sequence) : 
FRS : full resolution sequence 
WD : wavelet decomposition 
LRS : low resolution sequence 

MC-3DSA : motion-compensated 3D subband analysis 
30 LRD : low resolution decomposition 

HS : high subbands 

U-HFSS : union of the three high frequency spatial subbands of a frame 
FR-3D-SPIHT : full resolution 3D SPIHT 
OCB : output coded bitstream. 
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The corresponding decoding scheme, depicted in Fig.4, is symmetric to this encoder (in 
Fig.4, the additional references are the following : 

MC-3DSS : motion compensated 3D subband synthesis 

HSS : high subbands separation 

FRR : full resolution reconstruction). 
To enable spatial scalability, the high frequency spatial subbands just have to be cut as in the 
usual version of the 3DS codec, the decoding scheme of Fig.4 showing how to naturally 
obtain the low resolution sequence. 

Then, for coding the high spatial subbands, two main solutions are proposed, 
the first one without MC, and the second one with MC. 



A) Without MC 

In the first solution, the high subbands simply correspond to the high 
frequency spatial subbands of the original (full resolution) frames of the GOF in the wavelet 
decomposition. Those subbands allow the reconstruction at full resolution at the decoder. 
Indeed, the frames can be decoded at the low resolution. However, these frames correspond 
to the low spatial subband in the wavelet analysis of the original frames. Hence one has 
merely to put the low resolution frames and the corresponding high subbands together and 
apply a wavelet synthesis to obtain' the full resolution frames. But now, where and how to put 
those high subbands in order to optimize the 3D-SPIHT encoder ? In a MC scheme for a 3D 
subband encoder, the low temporal subbands always look like one of the original frames of 
the GOF. As a matter of fact : 

L=-~ , [A + MC(B)] . (3) 

so L looks like A. Consequently, the high spatial subband of A should be placed with the low 
resolution decomposition corresponding to L. This approach (reordering of the high spatial 
subband in the case of forward motion compensations) is illustrated in Fig.5, where DWTh 
denotes the high frequency wavelet filter and the coefficients Cj t are multiplication 
coefficients. The way to define Cj t is described later. 

However, the motion compensation in the 3D subband structure can be either 
forward or backward (it has even been shown that alternate directions improve coding 
efficiency. The following algorithm, in which the notations are : 

. jt : temporal decomposition level (0 for the full frame-rate, 
jtjnax for the lowest frame-rate) 
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. t : 0 for the low temporal subband, 1 for the high one 
. nf : subband index at temporal level jt 
. mejiir_desc_Jree : a byte that describes the ME directions 
used at a given temporal level jt (the LSB describes 
the direction of the first ME/MC, 0 means "forward", 
1 means "backward"), 

makes the link between a frame GOF_index in the GOF and the spatio-temporal subband 
{jt;n;t} which resembles it most, depending on the Motion Estimation Direction Description 
Tree. 
UTntS 

STlocationToGofIndex(MEDixectionDescriptionTree me__dir__desc_tree, UInt8 
jt_max,UM8jt,UIht8 nf,UInt8t) 
{ 

UlntS gof_index=0 ; 
UInt8 direction ; 
UInt8 j,n_sb ; 
UInt8 sign ; 

gofjndex - nf«jt ; 

sign = 1 ; 
n_sb = nf ; 

for(j^t-l;j>=0;j~) 
{ 

direction = l«n_sb ; 

if(t=0) 

sign=0 ; 

direction &= me_dir_desc_tree.aui8_level[j] ; 
direction »= n_sb ; 
if (sign) 
{ 

direction = ! direction ; 
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} 

n_sb = (n_sb«l) + direction ; 



direction «=j ; 
gof_index = direction; 

} 

return(gof_index) ; 



} 



The way to define the coefficients cj t is now described (in Haar filter case). Let 
a be the coefficient used in the temporal 2-tap Haar filter. In the conventional 3D subband 
scheme, one has : 

fL = a*(A+MC l (B) 
[H = a*(MC(A)-B) 

If, in the present scheme, one uses c Jt = a Jt for the high spatial subbands, then 
it is still meaningful to use temporal scalability. Indeed : 

DWT L (L) = a * ( DWT L (A) + MC" 1 (DWT L (B)) ) 
DWT H (L) = c jt *(DWT H (A)) 

= a * ( DWT H {A + UpSampl^UC' 1 (DWT L (B))] ) 
and : 

fDWT L (H) = a*( DWT L (B)-MCPWT L (A)] ) 
\DWT h (H) = a * ( DWT H (B) ) 

where UpSample refers to the picture upsizing using wavelet filters. For the reconstruction at 
a lower frame rate, only the low temporal subband is synthesized : 

L=^i-DWT" I [DWT(L)] 

= ^*( A+UpSampletMC^CDWTLCB))] ) 



Finally, the reconstructed frames at each temporal level will tend to look like a motion- 
compensated average of the "reference" original frame and a blurred version of the other one 
(up-sampled version of the downsized frame), whereas in the current version of the 3D 
subband codec this blur is not introduced. Improving spatial scalability at the expense of 
adding blur in the temporal scalability is however a worthy step. 
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B) With MC 

As using MC in every subband does not allow a reconstruction with no drift, it 
is possible, as depicted in Fig.6, to partially use MC to construct the high spatial subbands 
(which is better in tenns of coding efficiency) and still be able to reconstruct every resolution 
5 (in Fig.6, the additional references are the following : 

ME/MC : motion estimation/motion compensation 

PRE : prediction error). 
Instead of directly using the high frequency spatial subbands of the wavelet decomposition, a 
wavelet decomposition is carried out on a prediction error obtained from the MC performed 
10 on the full resolution sequence and reusing for instance the motion vectors of the low 
resolution. 

The solution is to define : 

|DWT h (L) = c ^ * (DWT H (A)) 

[D WT H (H) = c j( * D WT H (B - MC(A)) 

It can be noticed that the MC is only used in the high temporal subband : A is first 

1 5 reconstructed at the full resolution thanks to the low temporal subband, and then used to get 
frame B with MC thanks to H. The coefficients Cj t are chosen as previously. Said MC at full 
resolution can be performed either by merely upsampling the low resolution motion vectors 
(which has the advantage of introducing no other motion vector overhead) or by refining 
these upsampled low resolution vectors (which costs some additional transmission bits but is 

20 more efficient in terms of texture coding). 

It must be understood that the present invention is not limited to the 
aforementioned embodiments, and variations and modifications may be made without 
departing from the spirit and scope of the invention. There are numerous ways of 
implementing functions of the method according to the invention by means of items of 

25 hardware or software, or both, provided that a single item of hardware or software can carries 
out several functions. It does not exclude that an assembly of items of hardware or software 
or both carry out a function, thus forming a single function without modifying the method in 
accordance with the invention. Said hardware or software items can be implemented in 
several manners, such as by means of wired electronic circuits or by means of an integrated 

30 circuit that is suitable programmed. The integrated circuit can be contained in a computer or 
in an encoder or decoder and comprise a set of instructions, contained, for example, in a 
computer programming memory or in an encoder or decoder memory and causing the 
computer or the decoder to carry out the different steps of the methods according to the 
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invention. This set of instructions may be loaded into the programming memory by reading a 
data carrier such as, for example, a disk. A service provider can also make the set of 
instructions available via a communication network such as, for example, the Internet. 



